Fine-tuning the accretion disk clock in Hercules X-l 
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ABSTRACT 

RXTE ASM count rates from the X-ray pulsar Her X-l began falling consistently during the 
late months of 2003. The source is undergoing another state transition similar to the anomalous 
low state of 1999. This new event has triggered observations from both space- and ground-based 
observatories. In order to aid data interpretation and telescope scheduling, and to facilitate 
the phase-connection of cycles before and after the state transition, we have re-calculated the 
precession ephemeris using cycles over the last 3.5 years. We report that the source has displayed 
a different precession period since the last anomalous event. Additional archival data from 
CGRO suggests that each low state is accompanied by a change in precession period and that the 
subsequent period is correlated with accretion flux. Consequently our analysis reveals long-term 
accretion disk behaviour which is predicted by theoretical models of radiation-driven warping. 

Subject headings: accretion, accretion disks — instabilities — ephemerides — (stars:) binaries: close — 
X-rays: binaries — X-rays: individual (Her X-l) 



1. Introduction 

Her X-l is an eclipsing X-ray pulsar that dis- 
plays variability on spin (1.24s), orbital (l-7d) and 
super-orbital (35d) timescales (Tananbaum et al. 
1972). The companion star is of early-F or late-A 
type and fills its Roche lobe, resulting in accretion 
mainly by Roche lobe overflow through the inner 
Lagrangian (Li) point (Leahy & Scott 1998). 

The 35d period is revealed in X-rays as an ab- 
sorbing column of variable depth passing in front 
of the accretion source (Scott, Leahy & Wilson 
2000), In the UV and optical, an intrinsic flux 
variation is observed on the same period (Gerend 
& Boynton 1976; Leahy & Marshall 1999). The 
long-standing model is of a warped accretion disk, 
viewed close to edge-on (i ~ 90 deg) precessing 
in a retrograde direction around the neutron star, 
(Katz 1973; Roberts 1974). The disk atmosphere 
provides the column and casts an X-ray shadow 
over the inner face of the companion star, causing 
optical variability on 35d timescales. 
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In X-rays, the 35d cycle is characterized by a 
rapid turn-on, a "main-on" state which varies in 
intensity greatly from cycle to cycle, and then de- 
cays back to an "off" state at ~ 1% of the peak 
main-on flux. This is followed by a rise to a weaker 
"short-on" state and then another low (Tanan- 
baum et al. 1972; Gorecki et al. 1982). Superim- 
posed on this cycle are eclipses of the neutron star 
by the companion once per orbit and dips that oc- 
cur almost on the beat period between orbit and 
precession (Giacconi et al. 1973), which are be- 
lieved to be due to localized structure in the disk, 
generated by the impact of the ballistic gas stream 
from the L\ point (Crosa & Bonton 1980; Schandl 
1996). The turn-on occurs at only two orbital 
phases, </> or b — 0.2 or or b — 0.7, apparently at 
random. Consequently it has often been claimed 
that precession is not strictly periodic (Ogelman 
1987), and precession cycles generally have dura- 
tions of either 20, 20.5 or 21 orbital cycles, lead- 
ing some authors to suggest that the precession 
and orbital cycles are related physically (Scott & 
Leahy 1999). 

To date, it has not been possible to determine 
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Fig. 1. — Top panel: Her X-l CGRO BATSE (20-50 keV) one-day average count rates from Earth occulta- 
tions between 1991-1999. Middle panel: RXTE ASM (1-12 keV) one-day average count rates from Her X-l 
from 1996-2004. Bottom panel: ASM Hardness ratio averaged over each 35 day cycle. Band A = 1.3-3 keV, 
band B = 3-5 keV and band C = 5-12 keV. Horizontal lines represent count rate thresholds adopted in the 
— C analysis. Dwells over the two thresholds are plotted on the O — C diagram of Fig. 2 
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the engine behind disk precession observationally. 
From a theoretical perspective, the companion 
star can force a tilted disk to precess in the outer 
regions closest to the donor star (Larwood 1998). 
However tidal forces do not account for the mech- 
anism that tilts the disk. Mechanisms that can si- 
multaneously warp structure and drive precession 
are wind (Schandl & Meyer 1994) and radiation 
(Pringle 1996; Wijers & Pringle 1999) pressure, 
pushing down on the disk plane. The dominant 
engine behind this pressure is accretion flux from 
the neutron star which irradiates the disk atmo- 
sphere externally. While the simplifying assump- 
tions used in these theoretical calculations prevent 
a detailed comparison with observation, the suc- 
cess of the models in producing warped disks with 
the observed precession rate indicates that some 
combination of tides, wind and radiation pressure 
are capable of driving the Her X-l clock. 

Her X-l has occasionally displayed consecutive 
main-on states of significantly lower flux than av- 
erage, e.g. 1993 (Vrtilek ct al. 1994). Also, be- 
tween 1999-2000, the source missed a large num- 
ber of consecutive cycles altogether (Parmar et al. 
1999; Coburn et al. 2000; Vrtilek et al. 2001; Still 
et al. 2001). See also Parmar et al. (1985) for the 
description of a 1983 event prior to the epoch of 
all-sky monitoring surveys. During these times, 
the optical flux does not dip appreciably (Vrtilek 
et al. 2001). In addition to the decrease in main-on 
amplitude, the neutron star experiences an abrupt 
transition from the nominal state of pulse period 
spin-up to episodes of spin-down (Parmar ct al. 
1999). Although these events are not homoge- 
neous in terms of either duration or magnitude, 
they have been coined collectively the 'Anomalous 
Low States (ALS)'. Quite clearly the ALS are sig- 
nificant disk events; they point to the accretion 
disk changing state for a limited period of time 
before returning seemingly to the same preferred 
period and warp shape, albeit with an offset of 
several orbital cycles in the turn-on clock, accord- 
ing to Oosterbroek et al. (2001) and Manchanda 
(2003). In this paper, we show this picture to be 
incorrect. 

Rossi X-ray Timing Explorer (RXTE) All-Sky 
Monitor (ASM) count rates from Her X-l in the 
main-on state began falling consistently during the 
late months of 2003. Her X-l appears to be un- 
dergoing another state transition similar to the 



ALS of 1999-2000. Since a glitch in the preces- 
sion clock was reported after the last ALS (Ooster- 
broek et al. 2001), it is now prudent to recalculate 
the precession ephemeris using the last 3.5 years 
of data. Rather than a glitch, we report that the 
source has displayed a different precession period, 
P35, since the last ALS. Additionally, archival data 
from the BATSE experiment that was onboard the 
Compton Gamma Ray Observatory (CGRO) sug- 
gests that each ALS is accompanied by a change 
in P 35 . Contrary to previous reports (Oosterbroek 
ct al. 2001), these results indicate that the accre- 
tion disk does not return to the same state after 
each ALS. 

2. Results 

Fig. 1 plots the total BATSE and ASM count 
rate as averages of all dwells taken within a 24h pe- 
riod, weighted by uncertainty. Three ASM energy 
channels allow the source's spectral evolution to 
be monitored coarsely as a hardness ratio. The 
main-ons of the 35d cycle reveal themselves as 
maxima in the count rate time series, where each 
cycle peaks at 3-9 ASM count s _1 (1 Crab = 75 
count s _1 ). The stand-out features of this time 
series are the extended 1999 ALS between MJD 
51 195-51892 (Parmar et al. 1999), and the de- 
cline into a putative new ALS starting MJD 52 866 
(Boyd & Still 2004). The 1993 ALS reported by 
Vrtilek ct al. (1994) can be seen as two consecu- 
tive main-ons with faint maxima of < 0.05 BATSE 
count s" 1 between MJD 49 213-49282. 

The extended ALS also reveals itself in the 
hardness data. A soft spectrum in the ALS is con- 
sistent with that seen during normal low states of 
the 35d cycle (Scott & Leahy 1999; Clarkson et al. 
2003). There is some evidence that the hardness 
ratio returns to a normal level before the main-on 
cycle returns in count rate, and of long-term hard- 
ness evolution over the duration of ASM monitor- 
ing. For an analysis of the structure of ASM count 
rates, folded over the 35d cycle, see Shakura et al. 
(1998), Scott & Leahy (1999) and Leahy (2002). 
A more detailed analysis of 35d spectral evolution 
with GINGA is presented by Leahy (2001). 

A period search over the ASM data obtained 
after the 1999 ALS between MJD 51892-52866 
using the Lomb-Scargle method (Scarglc 1982) re- 
veals a best period of 35.10 ± O.Old, with a false- 
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Fig. 2. — Observed — Computed time series of the main-on states in the Her X-l precession cycle. Each point 
is a daily dwell in the time series that contains count rates above a threshold of 0.05 count s _1 (BATSE) or 
5 count s _1 (ASM). The reference ephcmcris adopted is T35 = 50986.7 + 34.79 E35. Dashed lines represent 
linear least-square fits to the O — C residuals between each ALS. Dotted vertical lines represent the start 
and stop times of each ALS. 



Table 1: Precession periods, P35, period derivative, P35, and maximum averaged main-on state count rates, 
Pnax, from four epochs of the Her X-l precession cycle. Cycle number is relative to MJD 50 986.7, for 
consistency with Fig. 2. 



Epoch 


Cycle 


MJD 


P35 


P S 5 


^max 




Start 


Stop 


Start 


Stop 


(d) 


(dd- 1 ) 


(count 




1 


-75 


-52 


48 377 


49 213 


34.637 + 0.010 


-(3 + l)xl0~ Y 


0.046 ± 0.002 a 




2 


-49 


-8 


49 282 


50 708 


34.977 ± 0.003 


(3±2)xl0~ 8 


0.064 ±0.001 a 


, 4.79±0.05 b 


3 


-5 


+6 


50 813 


51 195 


34.745 ± 0.021 


-(7 + 7)xl0~ 10 


0.042 ± 0.003 a 


, 3.66±0.07 b 


4 


+26 


+54 


51892 


52 866 


35.137 + 0.007 


(1 ± l)xl0" 7 




5.00±0.06 b 



a CGRO BATSE 
b RXTE ASM 



4 



alarm probability < 10 6 . This compares directly 
with an identical search over the pre-1999 ASM 
data, MJD 50 813-51 195, of 34.79±0.03d. Clearly 
there is dichotomy here. To qualify this, we con- 
struct an "Observed — Computed" (O — C) curve 
of precession cycles (Fig. 2). Since turn-on times 
are not strictly periodic, the simplest reference 
point for this analysis is the mid-point of each 
main-on state, T ma ; n . Any data point > an ar- 
bitrary threshold of 0.05 (BATSE) and 5 (ASM) 
count s~ x is plotted on the O—C diagram. So each 
individual main-on is represented by several data 
points in the diagram. Scatter is a relative mea- 
sure of the duration of main-on states, although 
ASM epochs cannot be compared directly with 
BATSE epochs in this way because of the differ- 
ent flux thresholds. Our reference cphcmcris, mea- 
sured from all ASM data before the 1999 ALS, is 
T main = MJD 50986.7 + 34.79 £35, where E 35 is 
an integer. The quantity subtracted from BATSE 
and ASM times to produce the O — C time series 
is the nearest T ma i n to each daily dwell. 

There are at least four epochs of different pre- 
cession period evident in Fig 2. There is no evi- 
dence for episodes of high P35 , except at discrete 
times that correspond to each ALS. The one ex- 
ception is at MJD 50 750. There have been no re- 
ports of an ALS previously for this date, however 
we suggest that the low count rates and period 
change at this epoch qualify it as such. There 
is an unfortunate gap in the published Her X-l 
pulse period history (Parmar et al. 1999) during 
this epoch where one would expect an episode 
of spin-down from the neutron star. Individual 
Lomb-Scargle searches and quadratic fits to the 
O — C residuals provide independent, and con- 
sistent, measures of P35 for each epoch. O — C 
fits yield the smallest statistical error and we pro- 
vide these in Table 1. Measured values of P35 are 
small during all epochs and these are also included 
in Table 1. Best fit ephemerides are plotted over 
the O — C diagram as dashed lines, although we 
have made no attempt to align the ephemerides at 
the turning points. Based on our arbitrary count 
rate threshold, there are other candidates for ALS, 
particularly during the early BATSE era at MJD 
48 600. However, there is no evidence within the 
O — C residuals or pulse history diagrams to con- 
firm this. 

Epoch- folding of the data over P 35 , shows that 



turn-ons typically occur 0.47T in phase before the 
mid-point of main-on, as defined by our thresh- 
old criteria, although there is an uncertainly of 
a few days because the turn-ons and maxima are 
decoupled. As discussed in Sec. 1, turn-ons ap- 
pear not to be defined by precession period alone, 
also being dependent on orbital phase. Therefore, 
while we provide a precession ephemeris below for 
the turn-on times of epoch 4, T on , we must allow 
a systematic uncertainty of ~ 1 orbital cycle (~ 
2d): 

T on = MJD 51 821(2) + 35.10(1)P 3 5 (1) 

Fig. 3 plots the time series, separated by epoch 
as defined in Table 1, and folded over the rele- 
vant precession period. The data are averaged into 
50 phase bins, where ^35 is the precession phase. 
^35 = is defined as the turn-on time. There 
are variations in the main-on structure from epoch 
to epoch, most noticeably in the peak flux of the 
main on state. 

In Fig. 4 we plot main-on count rates, F max , 
versus precession period, P35. P max is measured 
as the peak of a Gaussian fit of the folded light 
curves in Fig. 3 between 35 = 0.0-0.4. We find a 
positive correlation between P 35 and P max in both 
the BATSE and ASM data, although the sample 
is, of course, very small, and the different energy 
bands of the two instruments means that they can- 
not be compared directly. 

3. Discussion 

Clarkson et al. (2003) noticed the period vari- 
ation in the Her X-l precession cycle at MJD 
50 750, but interpreted it as a phase shift, suggest- 
ing that this event resulted eventually in the 1999 
ALS, 400d later. From a longer baseline of data, 
we offer the alternative interpretation that each 
ALS is accompanied by a "simultaneous" , change 
in P35. The timcscalc of state changes must be 
just a few precession cycles. These are followed by 
epochs where P35~ 0, indicating long-lived, sta- 
ble disk structure. Extrapolating the quadratic 
O — C fit for epoch 4, we find that an instanta- 
neous P35 would have occurred towards the begin- 
ning of the 1999 ALS. A Gabor transform analysis 
(i.e. a short-time Fourier transform with a sliding 
window; Hcil & Walnut (1989), appropriate for 
uncovering evolution in frequency and amplitude, 
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35d phase (0 35 ) 

Fig. 3. — One-day BATSE and ASM dwells from each of the four epochs, folded over the epoch-dependent 
precession period and averaged into 50 phase bins. The epoch ranges are defined in Table 1. The 35d cycle 
is plotted twice for the interest of clarity. 
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Fig. 4. — The time-averaged main-on BATSE and ASM count rates between cf> 35 = 0.1-0.4 as a function 
of precession period, P35. Since the two instruments are sensitive to different energy ranges we have not 
renormalized one to the other, but have multiplied the BATSE count rates by 10 2 to incorporate them on 
the same plot. No attempt is made to correlate BATSE fluxes with ASM fluxes due to the discrepent energy 
ranges and instrument characteristics between the two detectors. 



resolves the ASM light curve of Her X-l into a 
series of distinct episodes of constant period, con- 
sistent with the measurements reported in Table 1. 
However, due to the smaller sampling, the period 
uncertainties resulting from the Gabor transform 
analysis are much larger than those reported in Ta- 
ble 1. This explains why the three ASM precession 
epochs were not distinguishable in the dynamical 
power analysis of (Clarkson et al. 2003). 

From Fig. 4, we can speculate that the main- 
on flux and precession period are related. At 
the peak of the main-on state, the absorbing col- 
umn in front of the pulsar is consistent with the 
galactic column (dal Fiume et al. 1998), which 
does not contribute to the BATSE or ASM energy 
bands. It is reasonable, therefore, to assume that 
-Fmax measures the accretion flux from the pulsar 
(Scott, Leahy & Wilson 2000). A relationship be- 
tween P35 and -F m ax is a natural prediction from 
models of radiation- and wind-driven precession 
(Schandl & Meyer 1994; Wijers & Pringle 1999), 
however an attempt to correlate the two quantita- 
tively would require radiative and hydrodynamic 



models of enormous complexity. Wijers & Pringle 
(1999) show that the precession timescale (Pa) is 
governed by: 



where Rd is a characteristic accretion disk radius, 
Sd a surface density and is the strength of the 
radiation field from the disk. If the radiation field 
from the disk is dominated by irradiation from 
a central source then P max oc Since we find 
Pmax oc P35, Eqn. 2 is not consistent with observa- 
tion. However disk warp is self-limiting; brighter 
than some threshold L<j,max, the warp becomes so 
tightly wound that the disk is self-occulting and 
radiation pressure is reduced. In this phase space 
the precession period will be directly-correlated 
with the neutron star flux, because the surface 
area on the disk incident to the pulsar decreases 
with increasing warp. The simplest picture for 
Her X-l, therefore, is that the pulsar luminosity 
is close to the stability limit. This is consistent 
with the predictions from analytical modeling by 
Ogilvie & Dubus (2001). 
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A corollary from Eqn. 2 is that it is of no sur- 
prise that the precession period varies from epoch- 
to-epoch. In fact, it is predicted directly from 
theory (Wijers & Pringlc 1999; Ogilvie & Dubus 
2001). If structure extends to the outer disk edge, 
any alteration in disk tilt will result in the deposi- 
tion of mass by the ballistic stream at a different 
range of accretion disk radii (Lubow 1989; Schandl 
1996). The resulting redistribution of surface den- 
sity over the disk will yield a new equilibrium con- 
figuration for the new epoch, according to Eqn. 2. 

We have determined P35 from the entire time 
series rather than the interval between consecu- 
tive main-on turn-ons. Since turn-ons occur at 
only orbital phases <j) ib— 0.2 and 0.7, it is nat- 
ural to assume that structure locked in the rota- 
tional frame of the binary is partly responsible for 
turn-on times (Schandl 1996). Measuring inter- 
vals between consecutive turn-ons would provide 
the same results as Sec. 2 only if a large number of 
cycles were sampled. As we have shown, P 35 has 
not remained stable for more than 40 consecutive 
cycles in the BATSE and ASM era. Consequently 
measuring turn-on intervals provides a more un- 
certain measure of P35. One interesting direction 
for future work would be to determine whether 
turn-on times can be predicted from P 35 , i.e., does 
the time of turn-on correspond to whichever or- 
bital phase, 0.2 or 0.7, immediately follows an un- 
determined reference phase within the stable pre- 
cession cycle? Of course, the reference will change 
after each ALS. 

4. Conclusions 

We have redefined the accretion disk precession 
ephemeris for the X-ray pulsar Her X-l between 
the end of the 1999-2000 ALS and the start of 
the 2003-2004 ALS event. Relative to the most 
recently published work on cycle timing (Man- 
chanda 2003), the new ephemeris provides a cor- 
rection to the turn-on time of the precession cycle 
of ~ 14d at the time of writing, critical to research 
teams and observatory schedulers during the cur- 
rent ALS. A period analysis of the CGRO BATSE 
and RXTE ASM time series since 1991 reveals four 
epochs of stable precession period. The preces- 
sion epochs are each separated by an ALS and 
the change of period appears to be instantaneous 
within the resolution of the O — C diagram, i.e. 



a few precession cycles. Testable predictions are 
that Her X-l will return from the 2003-2004 ALS 
with a new precession period and epoch-averaged 
peak main-on flux. 

Results provided by the RXTE ASM teams at 
MIT and at the RXTE SOF and GOF at NASA's 
GSFC, and the CGRO BATSE instrument team 
at NASA's MSFC. 
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